Graphene superlattice with periodically modulated Dirac gap 
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Graphene-based superlattice (SL) formed by a periodic gap modulation is studied theoretically 
using a Dirac-type Hamiltonian. Analyzing the dispersion relation we have found that new Dirac 
points arise in the electronic spectrum under certain conditions. As a result, the gap between 
conduction and valence minibands disappears. The expressions for the position of these Dirac 
points in k-space and threshold value of the potential for their emergence were obtained. At some 
parameters of the system, we have revealed interface states which form the top of the valence 
miniband. 

PACS numbers: 71.10.Pm, 73.21.-b, 81.05.U- 



I. INTRODUCTION 

During for the last years extremely much attention was 
paid to the electronic properties of graphene (see Ref.l for 
the review). Such interest results, in particular, from the 
fact that physics of the low-energy carriers in graphene is 
governed by a Dirac-type Hamiltonian. The band struc- 
ture of an ideal graphene sheet has no energy gap. As 
a consequence, Dirac electrons become massless, and re- 
veal unusual properties such as perfect transmission at 
normal incidence through any potential barrier(the Klein 
paradox^), trembling motion (or Zitterbewegungi), etc. 
Because of the Klein effect, an electrostatic potential can- 
not confine electrons in graphene. This makes difficult 
to use graphene in electronic devices £ To overcome this 
problem, one needs to create a gap in graphene. An- 
other possible way for confining Dirac fermions is the use 
of spatially inhomogeneous magnetic field £ The gap can 
be induced by substrate or strain engineering as well as 
by deposition or adsorption of molecules on a graphene 
layer. For example, two carbon sublattices of graphene 
placed on top of hexagonal boron nitride (/i-BN) become 
nonequivalent due to their interaction with the substrate. 
The band-structure calculations within the local-density 
approximation for this system gives a gap not less than 
53 meVj£ A hydrogenated sheet of graphene (graphane) 
is a semiconductor with a gap of the order of a few eVfi 

Besides, it is possible to modulate spatially the gap 
(i.e. the particle's mass) in graphene. It was shown that 
the spatial mass dependence leads to suppression of Klein 
tunneling and induces confined states The required 
gap modulation can be created, for instance, in graphene 
placed on a substrate fabricated from different dielectrics. 
It is also possible to use for this purpose an inhomoge- 
neously hydrogenated graphene or graphene sheet with 
nonuniformly deposited Cr03 molecules. 

Correspondingly, one can fabricate different graphene 
heterostructures with the gap discontinuity. In particu- 
lar, graphcnc-bascd superlattice (SL) can be formed due 
to the periodic modulation of the band gap. The model 
of such a superlattice has been proposed and discussed 
earlier 



In our work we analyze the electronic states in 
graphene-based SL with periodically modulated gap and 
relative band shift (potential). We consider the simplest 
case where the gap and potential are piecewise constant 
functions, and compute the SL band structure at various 
relationships between the gap and band shift. We have 
found in the SL spectrum new Dirac points whose po- 
sitions in k-space were determined. We also show that 
interface states can arise in the gap-induced SL. For this 
case the relevant conditions and system parameters were 
obtained. 



II. DIRAC PARTICLE IN GRAPHENE 
SUPERLATTICE 

Let us consider one-dimensional (ID) graphene super- 
lattice with period / formed by position-dependent gap 
and band shift. As was shown in Ref.ll, such structure 
can be realized, e.g., on the base of graphene deposited on 
a strip substrate combined from silicon oxide and /i-BN 
(Fig.l). The SL electronic structure in the vicinity of K- 
point of the Brillouin zone is described by the Dirac-like 
Hamiltonian 



H = v F pa + V(x)l + A(x)a z 



(1) 



where p is the momentum operator, cr; are the Pauli ma- 
trices, 1 is a unit 2x2 matrix, vp « 10 s cm/s is the 
Fermi velocity, and A(x), V(x) are periodic functions 
equal to A and V, respectively, at a < x < I, and zero 
at < x < a. Here, the potential V defines the shift of 
the forbidden band center in the gapped graphene with 
respect to the Dirac point in the gapless graphen o 11 ' 12 
(see Fig.l). 

The Dirac equation 



HV(x,y) = EV(x,y) 



(2) 



admits the solutions of the form ^(x, y) — 
exp(ik y y)ty (x) , where the two-component spinor 
envelope function ^>(x) satisfies the equation 



i— = h(x)^(x) 
ax 



(3) 
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FIG. 1: (Color online) Top: graphene layer on the striped 
substrate composed of silicon oxide and hexagonal boron ni- 
tride. Bottom: schematic diagram showing the electronic en- 
ergy spectrum in graphene SL. 



with 



V(x)-E-A(x) 

The formal solution of Eq.(3) is 

^(x) = ?k x exp ( -i [ h(x 1 )dx 1 )^(x ), (5) 



-ik 



•i 



where ik x is the spatial ordering operator ! 13 ' 14 This ex- 
pression can be simplified as 



*(.t) = exp(-i(ac - x a )h)^(x ) 



(6) 



if both points x and xq belong to the space-homogeneous 
region. 

It is convenient to define the matrix 



t(x — xq) — exp(— i(x — xo)h). 



(7) 



Then = t(l - o)t(a)*(0) = T*(0). Performing the 
ordinary expansion of exponential function in power se- 
ries and using Eq.(4) one gets: 

i(x - Xq) = 

(k„ ■ ■ Vtx)-E-A(x) \ 

cosa+j^sma -isma \ j VfK(x) ' \ 

■ ■ V(x)-E+A(x) k v ■ ' \°) 

- tSmQ hv F K(x) COSa -#) SmQ / 

where K{x) = ^J{{V{x) - E) 2 - A 2 (x)) / {hv F f - k 2 yl 

a = (x — xo)K(x). Note, that det(f) = 1 and, conse- 
quently det(T) = 1 as well. This equality and the Bloch 
condition = ^(0) exp(ifeZ) (here, k is the Bloch wave 
vector) yield the dispersion relation 2cos(fcZ) = Tr(T) 
for the ID graphene-based SL. Using Eq.(8) one can find 
Tr(T) = Tr(t(l - a)t(a)). Accordingly, the resulting dis- 
persion equation reads 

cos(kl) — cos(k x a) cos (q x (l — a)) + 
EV-{hv F k x ) 2 , 



with k x = K(+0), q x = K(l - 0). At (V - E) 2 - A 2 < 
(TiVFky) 2 , the wave number q x is imaginary and the al- 
lowed energies are given by Eq.(9) with q x replaced by 

Eq.(9) for the allowed energies was obtained also in 
Refill via wave function matching. However some results 
of this work concerning the miniband structure and exis- 
tence of the interface states seem questionable. Below we 
discuss these issues in detail. Note that at A = 0, Eq.(9) 
coincides with the one found for single-layer graphene in 
a periodic piecewise constant potential V(x)M ~— 



III. ELECTRONIC STRUCTURE 

The SL band structure depends on the system parame- 
ters as well as on the y-component of the wave vector. As 
seen from the Eq.(9), the dispersion relation is invariant 
with respect to the simultaneous replacements E — > —E, 

V — > — V. Therefore, in what follows we shall consider, 
for dcfinitcness, only nonnegative values of the relative 
band shift: V > 0. When V — 0, the energy spec- 
trum is completely symmetric related to the value E = 
corresponding to the original Dirac point in the gapless 
graphene. In this case the two first minibands symmet- 
rically situated above and below the original Dirac point 
are the conduction and valence ones, respectively. As V 
increases, the conduction and valence minibands gradu- 
ally shift up. Below we shall concentrate on these two 
minibands only, assuming the Fermi level to be in be- 
tween (or, at least, within one of the minibands) at any 
V. 

The electron and hole energies as functions of k at 
k y = 0, A = 26.5 meV and for different values of the 
potential V and for different widths a are plotted in Fig. 2. 
The parameters we choose are appropriate for graphene. 
Dots and dashed (red) line correspond to V = and 

V = 24.2 meV, respectively and a — 1/2. It is seen that 
when V becomes nonzero, the electron band turns out 
to be narrower than the hole one. The electron and hole 
energies for a = 1/6, V — 24.2 meV are shown with solid 
(blue) line. We can see that, as the gapped graphene 
fraction in the SL increases, the electron-hole minigap 
increases too. Nevertheless, in any case, at k y = the 
minigap cannot exceed the gap value 2A»i£ Indeed, it 
is clear that if the whole graphene layer is gapped, i.e. 
a = 0, the forbidden miniband should be 2A (at k y =0). 

Figure 3 illustrates the dependence of the electron and 
hole energies on k y at different V and A, and k = 0, a = 
1/2, 1 = 60 nm. We show only semiaxis k y > because of 



the symmetry k y 



-k y in the dispersion low (9). The 



(Tiv F ) 2 k x q x 



■ sin(fc x a) sin (q x (l — a)) , (9) 



energies in the figure are counted from the minigap cen- 
ter (whose position determining at k = k y = depends 
on V) for V = 24.2 meV (short-dashed line) and V = 110 
meV (long-dased line) at A = 26.5 meV. In other cases 
the energy origin coincides with the contact-point energy. 
Our calculations demonstrate a monotonous expansion of 
the forbidden miniband with \k y \ increasing (short- and 
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FIG. 2: (Color online) Two low-energy minibands of the SL 
spectrum with k y — 0, I = 60 nm, and A = 26.5 meV. The 
dotted (violet) line corresponds to the barrier height V = 
and a — 1/2. The solid (blue) and dashed (red) lines corre- 
spond to the width of gapless graphene a = 1/6 and a — 1/2 
respectively; V — 24.2 meV in both cases. 



long-dashed (blue) lines in Fig. 3) until V becomes greater 
than some threshold value V c (for chosen set of parame- 
ters V c = 143.26 meV). When V = V c (thick solid line in 
Fig. 3), the electron and hole energy branches touch each 
other at k = k y = closing the minigap. When V ex- 
ceeds V c , the minigap at k = k y = opens, but two extra 
Dirac points appear in symmetric positions on the fc y -axis 
(thin solid line in Fig. 3), and never then vanish. Thus, 
beginning with V = V c the gap-induced SL becomes gap- 
less. Further increase of V results in appearance of new 
additional Dirac points, which originate from k y = at 
certain values of the potential V (see below). 

Such a behavior is not unique. As was shown in Rcfs 
18-20 for gapless graphene, the new Dirac points can 
emerge in the presence of a sinusoidal or squarewell SL 
potential. A simple analytical expression for the location 
of these points in k-space has been obtained by different 
research groupsi 14 i 16 

It is clearly seen from Fig. 3 that extra Dirac point 
located at k y l ~ 6.3 for A = (dots) remains also in the 
case where the SL consists of both gapless and gapped 
graphene. However, the position of the Dirac point in this 
case slightly shifts towards zero. To find the location of 
the Dirac points in k-space we note that k x {E,k y ) and 



q x {E,k y ) coincide at 



E (V) 



V 2 - A 2 
2V 



(10) 



Correspondingly, Eq.(9) at a = 1/2, k = and E = E 
turns into 

2 kj EpV - Ej + h 2 v 2 F k 2 y , 2 k x l 
1 = C ° S - + (h^kj 2 Sm (U) 

Obviously this equation is satisfied if k x l = 27m (n^0), 
which leads to 



(Eol_\ 
\hv F ) 



- (27m) 5 



(12) 
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FIG. 3: (Color online) fcy-dependence of the electron energies 
in the conduction and valence minibands at k — 0, I — 60 nm, 
a = 1/2, A = 0, and V = 196 meV (dots (red)). For other 
curves A = 26.5 meV and V = 24.2 meV (short-dashed line); 
V = 110 meV (long-dased line); V = 143.26 meV (thick solid 
line), and V = 196 meV (thin solid line). The zeroth energy 
corresponds to the minigap center for the dashed lines, and 
to the contact points for all the others. 



where n is integer. 

One can check that for a — 1/2, the right side of Eq.(9) 
has local maximum (equal to unity) at E = Eq and k y de- 
fined by Eq.(12). This means coincidence of conduction 
and valence miniband edges. 

According to Eq.(12), in general case, when the ratio 
Eol/2Trhi>F is not integer, the number of Dirac points Nd 
symmetrically situated around k y = is given by 



N. 



D 



EqI 



2nhvf 



(13) 



where [...] denotes an integer part. One of the contact 
points is always located at k y = if EqI/%Vf = 2nj 
(j = 1,2,...). In this case the total number of Dirac 
points is odd: Nd = 2j — 1, and the dispersion along 
k y -axis is almost flat in the vicinity of the Dirac point 
situated at k y = 0: E — Eq ~ ±k 2 A/n 2 . Similar behavior 
of the dispersion curves takes place in SL based on a 
gapless graphene^ where, however, dispersion is defined 
by the third power of the y-component of the wave vector: 
E — Eq ~ ±k y A/n 2 . Note that, the coefficient before k 2 
in our case is proportional to A only if A <C 2ttTivf/1- In 
other case, its dependence on A becomes more complex. 

Thus, now we can write analytically the inequality for 
the system parameters corresponding to the formation of 
the gapless SL. Evidently, the gap separating the con- 
duction and valence minibands vanishes if 



(V 2 - A 2 )l 



> 1. 



(14) 



AttHvfV 

The case of the rigorous equality in Eq.(14) corresponds 



4 



> 




Since the left side of this expression is positive (see 
Eq.(15)) the allowed values of the energy should be nega- 
tive if V > and vice versa. It is not difficult to show that 
the inequalities (15) and (17) have the solutions when A 
and V differ from zero, A > |V|, and k 2 < k%, where 
fc 2 = A 2 {A 2 ~V 2 )/(hv F V) 2 . 

In Fig. 4 the hole miniband of the SL with a — 1/2, 
I = 60 nm is plotted for V = 11 meV, A = 22 meV, and 
kyl = 1. Upper narrow slice of the valence miniband de- 
picted by dots is formed by the interface states. Indeed, 
it is possible to verify that any energy value from this 
range obey the inequalities (15) and (17) for the above 
parameters. Dashed (red) line in the figure shows the 
interface energy level \E\ — hvp\k y \(l — V 2 / A 2 ) 1 / 2 for 
an isolated heterojunction.— 



FIG. 4: (Color online) Hole band of the spectrum of SL for 
a = 1/2, I — 60 nm. Top of the energy band depicted by the 
dots (red) corresponds to interface states. Dashed (red) line 
shows the interface energy level for an isolated heterojunction. 



V. CONCLUDING REMARKS 



to the threshold value of the potential: V — V c (the 
positive root of the quadratic equation). 

We do not discuss here the case a ^ 1/2. We ex- 
pect that Dirac points appear in this case too. However, 
the electron-hole miniband profile should be significantly 
asymmetric related to the Fermi level similarly to that 
takes place in SL based on a gapless graphene^ 



IV. INTERFACE STATES 

As was shown by Ratnikov and Silin 21 interface states 
(ISs) can exist in graphene-based hetero junctions. It 
was found that ISs result from the crossing of disper- 
sion curves of gapless and gapped graphene modifica- 
tions. However, the author of Ref.ll stated that appear- 
ance of IS minibands is impossible in a graphene-based 
superlattice. In contrast to this conclusion we found that 
the ISs can arise under certain conditions discussed be- 
low. Since the wave functions of these states behave as 
exponentials along x-asds, the wave vectors k x and q x 
should be imaginary, so that 

(hv F k y ) 2 >E 2 , (hv F ky) 2 + A 2 > {E -V) 2 . (15) 

In this case Eq.(9) transforms into 



cos(kl) — cosh(k' x a) cosh(q' x (l — a)) 



EV - E 2 + (hv F k y 



(hv F ) 2 k x q' x 



• sinh(k' x a) sinh(q x (l — a)), (16) 



with k' 



k% 



E 2 



Evidently the solution of this equation can exist only 



A 2 -(V-E) 2 



1.2 

K y 



if 



(HVFky 



E 1 < -EV. 



(17) 



We have considered the simple model of a one- 
dimensional SL in which the gap and potential profile 
are piecewise constant functions. In the framework of 
this model the dispersion relation for Dirac electrons was 
obtained, and the structure of low-energy minibands was 
investigated depending on the potential V and other pa- 
rameters of the SL. 

It was found that if V becomes greater than some criti- 
cal value V c , the electron- hole minigap turns out to be an 
essentially non-monotonous function of the y-component 
of the electron wave vector. In particular, when V > V c 
the forbidden miniband turns into zero at some values of 
k y , and new Dirac points arise symmetrically relatively to 
k y = in the electronic spectrum. As a result, the SL be- 
comes gapless. Appearance of new Dirac points at k y ^ 
is typical as well for SL induced by the potential mod- 
ulation of a gapless graphene i 14 i 16 i 18 ~ — Thus, one may 
conclude that such a reconstruction of the electron spec- 
trum in a certain extent represents a universal property 
of graphene-based superlattices. In the case where the 
widths of the gapless and gapped graphene strips in the 
SL are equal, we found the positions of the Dirac points 
in k-space and obtained the expression for the threshold 
value V c of the potential responsible for their appearance. 

Finally, we found that the interface states can exist in 
the gap-induced SL at certain conditions. 
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